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§ 1 .  Introduction . 


This  report  presents  an  analysis  of  the  structure  of  transient 
acoustic  waves  in  stratified  fluids  whose  densities  and  sound  speeds  are 
functions  of  a  single  depth  coordinate.  The  acoustic  field  in  such  a 
fluid  may  be  described  by  a  real-valued  function  u(t,x,y)  (the  acoustic 
potential  or  the  excess  pressure)  that  satisfies  the  wave  equation 

(1.1)  J^t  -  c2 (y)  p(y)  V  •  (p'VyJVu)  -  f(t,x,y) 

where  t  is  a  time  coordinate,  y  is  a  depth  coordinate,  x  *  (x1,x2)  are 
Cartesian  coordinates  in  a  horizontal  plane,  7  *  (3/3xt ,3/3x2,3/3y)  and 
f(t,x,y)  is  a  function  that  characterizes  the  wave  sources.  c(y)  and 
p(y)  are  the  variable  sound  speed  and  density,  respectively,  and 
p'1 (y)  -  l/p(y) . 

This  report  is  a  sequel  to  the  author's  report  [12]  on  "Spectral 
Analysis  of  Sound  Propagation  in  Stratified  Fluids."  That  work  contains 
a  spectral  analysis  of  the  acoustic  propagator 

(1.2)  Au  -  -cz(y)  p(y)  V  •  (p  ‘(y)7u) 

for  the  cases  of  an  unlimited  fluid  (-«  <  Xj,x2,y  <<*>),  a  semi-infinite 
layer  (-«  <  xltx2  <  «  and  0  <  y  <  “)  and  a  finite  layer  (-»  <  Xj,x2  <  <*> 
and  0  <  y  <  h  <  <*>) .  The  integration  of  (1.1)  below  is  based  on  the 
spectral  -analysis  of  [12]  and  the  notation  and  results  of  [12]  are  used 
throughout  this  report.  As  in  [12],  only  the  case  of  an  unlimited  fluid 
is  presented  in  detail.  The  modifications  required  in  the  second  and 
third  cases  are  described  in  §6  at  the  end  of  the  report. 
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The  analysis  of  the  structure  of  transient  acoustic  waves  in 
plane  stratified  media  was  initiated  by  the  author  in  1973  [7].  This 
work  dealt  with  the  special  case  of  the  Pekeris  profile  (propagation  in 
a  half-space,  p(y)  -  const.,  c(y)  piecewise  constant).  In  1974  the 
results  were  extended  to  the  symmetric  Epstein  profile  (p(y)  *  const., 
c  2(y)  ■  K  sech2  (y/H)  +  M)  [8].  The  general  Epstein  profile 
(p(y)  ■  const.,  c  “(y)  *  K  sech2  (y/H)  +  L  tanh  (y/H)  +  M)  was  treated 
in  1979  [11]  using  the  spectral  analysis  of  the  Epstein  operator  due  to 
Guillot  and  Wilcox  [2,  3].  A  preliminary  version  of  [11]  was  announced 
in  1978  [1,].  The  analysis  of  these  cases  was  based  on  explicit 
representations  of  their  normal  mode  functions  by  means  of  well-known 
special  functions.  Extension  of  the  analysis  to  larger  classes  of  strat¬ 
ified  fluids  had  to  await  the  extension  of  the  spectral  analysis  of  the 
acoustic  propagators  (1.2)  to  such  classes. 

The  class  of  stratified  fluids  treated  in  [12]  was  characterized 
by  the  properties 


(1.3) 


(1.4) 


(1.5) 


p(y)  and  c(y)  are  Lebesgue  measurable, 

0  <  Pm  <  P(y)  <  PM  <  0  <  <  c(y)  <cM<«, 


p(±«) | 


dy  <  », 


•  +oo 

|  c(y)  -  c  (±°°)  I  dy  <  », 

JO 


where  pm,  Pj^,  p(±°°),  c  ,  c^  and  c(±°°)  are  constants.  In  this  report  the 
results  of  [10]  and  [11]  on  transient  acoustic  waves  in  stratified  fluids 
are  extended  to  the  class  of  fluids  characterized  by  (1.3),  (1.4)  and 


► 
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(1.6) 


|p(y)  -  P(±°°)  |  <  C(±y) 


. 

|c(y)  -  c(±<=°) j  <  CCiy)-0 

j 


for  ±y  >  0, 


where  C  and  a  are  constants  and 


(1.7)  a  >  |  . 

It  is  clear  that  (1.6)  implies  (1.5).  It  will  be  seen  from  the  analysis 
below  that  (1.6)  could  be  replaced  by  other  order  conditions  at  y  ■  ±°°. 

It  is  not  known  whether  the  results  of  this  report  hold  for  the  entire 
class  of  luids  defined  by  (1.3),  (1.4),  (1.5). 

The  remainder  of  this  report  is  organized  as  follows.  In  §2  the 
normal  mode  expansions  for  unlimited  fluids  of  [12]  are  used  to  decompose 
transient  acoustic  fields  into  free  and  guided  components  and  to  obtain 
integral  representations  of  these  components.  The  behavior  for  large  t 
of  the  free  and  guided  components  is  derived  in  §3  and  §4,  respectively. 
In  §5  these  results  are  used  to  calculate  the  asymptotic  distribution  for 
large  t  of  the  wave  energy.  §6  presents  corresponding  results  for 
semi-infinite  and  finite  layers. 


§2.  Normal  Mode  Expansions  of  Transient  Acoustic  Fields. 

The  integration  of  the  wave  equation  (1.1)  will  be  based  on  the 
spectral  analysis  of  the  acoustic  propagator  A  in  the  Hilbert  space 
3C  ■  L2(R3,c  2(y)  p  1 (y)  dxdy) ,  as  developed  in  [12].  To  this  end  (1.1) 
is  interpreted  as  the  equation 

(2.1)  4tt  +  Au  =  f(t,*)»  t  e  R, 


for  a  function  t  •*  u(t,*)  €  3C.  The  wave  sources  will  be  assumed  to  act 
during  a  time  interval  [0,T],  so  that  supp  f  C  [0,T].  The  corresponding 
acoustic  wave  is  the  solution  of  (2.1)  that  satisfies  the  initial 
condition 

(2.2)  u(t,*)  ®  0  for  all  t  <  0. 

The  solution,  based  on  the  spectral  theorem  for  the  selfadjoint  operator 
A  >  0,  is  given  by  Duhamel's  integral 

(2.3)  u(t,*)  -  |  {A  1/2  sin  (t  -  x)Al/2}  f(T,*)dt,  t  >  0. 

*  o 

Indeed,  if  f  €  C([0,T],JC)  then  (2.3)  is  the  unique  "solution  with  finite 
energy"  of  [6],  while  if  f  6  C([0,T],D(A1/2  ))  then  (2.3)  is  the  "strict 
solution  with  finite  energy."  In  addition,  if  f  €  C( [0,T] ,D(A_1/2  ))  then 

(2.4)  u(t,x,y)  *  Re  (v(t,x,y)} 
where  v(t,*)  is  the  complex-valued  potential  defined  by 


6 


(2.5)  v(t,‘)  -  i  exp  {-it  Al/2}  A~1/2  j  exp  {ixA1/2}  f(t,*)dt. 


In  particular. 


(2.6) 


v(t,  •)  »  exp  {-it  A1^2}  h  for  all  t  >  T 


where 


(2.7) 


X 

i  A  ^2  |  exp  {it  A1/2}  f(t,*)  dt. 


The  initial  value  problem 


(2.8) 


+  Au 


0  for  t  >  0, 


(2.9) 


u(0) 


f  du(0) 
’  dt 


g 


can  be  treated  by  the  same  formalism.  Indeed,  if  f  6  D(Aly<2  )  and 
g  e  D(A-1^2  )  then  the  solution  of  (2.8),  (2.9)  is  given  by  (2.4),  (2.6) 
with  h  *  f  +  i  A_l/2  g  €  D(Al/2  );  cf.  [9,  Ch.  3]. 

The  integral 

(2.10)  E(u,K,  t)  -  |  {  | Vu  j  2  +  c”2(y)|§7j  }  p-1(y)  dxdy 

may  be  interpreted  as  the  energy  of  the  acoustic  field  u  in  the  set 
K  C  R3  at  time  t.  Moreover,  A  is  the  selfadjoint  operator  in  3C  asso¬ 
ciated  with  the  sesquilinear  form  A  on  JC  defined  by  D(A)  -  Lj(R3)  C  Jf 


and 
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(2.11) 


A(u,v) 


Vu  •  Vv  p  1 (y)  dxdy 

Jr3 


It  follows  from  Kato's  second  representation  theorem  [4,  p.  331]  that 
D (A 1/2  )  *  L2 (R3)  and  for  all  u  e  D(Al/z  ) 


(2.12)  I) A1^2  ufljf  *  A(u,u)  =  |Vu|2  p  1  (y)  dxdy. 

Jr3 

Hence  the  total  energy  satisfies 


(2.13) 


E (u,R3 , t)  =  BA172  u|  J 


9u 

at 


If  h  €  D(Al7z)  and  u  is  defined  by  (2.4),  (2.6)  then  a  simple  calculation 
shows  that 

(2.14)  E(u,R3, t)  ■  0A1/2  hS^.  for  all  t  >  T. 

The  analysis  of  the  structure  of  the  acoustic  potential  (2.4), 
(2.6)  presented  below  is  based  on  the  normal  mode  expansion  for  A  of 
[12,  §8].  The  orthogonal  projections  {P+,P_,P0,P1 ,P2, • • •}  in  3f  defined 
by  the  normal  modes  form  a  complete  family  that  reduces  A  [12, 
Corollaries  8.4  and  8.10].  Hence  the  same  is  true  of  the  orthogonal 
projections 


(2.15) 


Pf  -  P+  +  P_  +  P0 


and 


P 

g 


N0-l 

I  V 

k-1 


(2.16) 


This  provides  a  decomposition 


(2.17)  u(t,*)  =  uf(t,*)  +  ug(t,») 
into  orthogonal  partial  waves 

(2.18)  uf(t,‘)  =  Pf  u(t,-)»  ug(t,*)  =  Pg  u(t,*). 

The  function  u^,  called  the  free  component  of  u,  will  be  shown  to  behave 
for  large  times  like  a  diverging  spherical  wave.  The  function  u  ,  called 

O 

the  guided  component  of  u,  will  be  shown  to  be  approximately  localized 
near  one  or  more  of  the  planes  y  *  const,  where  c(y)  has  local  minima  and 
to  propagate  outward  in  horizontal  planes  like  a  diverging  cylindrical 
wave.  This  second  component  shows  the  profound  effect  of  acoustic  ducts 
on  transient  acoustic  waves.  It  is  absent  when  c(y)  has  no  local 


minima.. 


§3.  Transient  Free  Waves. 


The  normal  mode  expansions  of  [12]  are  used  in  this  section  to 
calculate  the  asymptotic  behavior  for  t  ■+•  “  of  the  free  component 
u^(tf*)  ■  u(t, •) .  The  principal  result  is  that  in  each  of  the  half¬ 
spaces  R_|_(d)  and  R*(d),  where  R2(d)  -  {(x,y)  :  ±(y  -  d)  >0},  uf(t,*)  is 
asymptotically  equal  to  a  wave  function  for  a  homogeneous  fluid  with 
propagation  speed  c(°°)  and  ci-00) ,  respectively.  It  is  this  behavior 
that  motivates  the  term  "free  component"  for  u^t,*)- 

It  will  be  assumed  that  the  total  acoustic  potential  u  satisfies 
u(t,*)  *  Re  {v(t,  •)}  where  v(t,*)  ■  exp  {-it  A^z}  h  and  he  D(A1/2) 

(see  §2).  The  corresponding  partial  waves  u^Ct,*)  *  u(t,*)  with  k  >  1 

satisfy  uk(t,*)  -  Re  (exp  (-it  A1/2)  P^h} .  This  follows  from  the  fact 
that  the  normal  mode  functions  (y , p)  are  real  for  k  >  1,  which  implies 
that  Pk(h)  »  Pk(h) .  It  follows  by  addition  that  ug(t,  •)  * 

Re  {exp  (-it  A1/2)  P  h}  *  Re  {v  (t,*)}  and  hence 
g  S 

[  uf(t,*)  -  Re  {vf(t,*)}. 


(3.1) 


vf(t,*)  -  exp  (-it  Al/2)  Pfh  -  exp  (-it  A1/2)  hf. 


The  normal  mode  representations 


(3.2)  vf(t,x,y)  -  <j)+(x,y,p,q)  exp  {-itX1/2  (p,q)}  h+(p,q)  dpdq 

JR3  ~ 

of  [12,  §8]  provide  the  starting  point  for  calculating  the  asymptotic 
behavior  of  u^t.x, y)  for  large  t.  The  integral  in  (3.2)  will,  in 
general,  converge  only  in  X.  For  brevity  the  Jf-lim  notation  of  [12]  is 
suppressed  in  this  report. 
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Equation  (3.2)  gives  two  representations  of  v^  corresponding  to 
the  two  families  <J>+  and  <£_.  The  calculations  below  are  based  on  the 
^-representation  which  has  been  found  to  yield  the  simplest  form  of 
the  asymptotic  wave  function.  It  will  be  convenient  to  introduce  the 
characteristic  functions  Xo  and  x_  of  the  cones  C+,  C0  and  C_  in 
(p»q)~space  [12,  (1.48)— (1.51) ]  and  to  decompose  h_  as 

(3.3)  h_(p,q)  =  £(p,q)  +  m(p,q)  +  n(p,q) 

/V  /V  A 

where  Z  ■  x+  h_,  m  *  Xo  h_  and  n  =  x_  h_.  The  corresponding  decomposi¬ 
tion  of  v ^  is 

(3.4)  vf  "  v£  +  vm  +  vn 
where 

* 

v^  -  exp  (-it  A1^2)  $>*& 

(3.5)  ■  v  *  exp  (-i  t  A1^2)  $*m 

m  — 

v^  *  exp  (-it  A1/2)  $*n 

The  behavior  for  t  **•  00  of  these  three  functions  will  be  analyzed 
separately. 

Behavior  of  v^.  The  partial  wave  v^  has  the  representation 

(3.6)  vJl(t,x,y)  -  Jc  <p_(x,y,p,q)  exp(-i  t  w+(p,q))  £(p,q)  dpdq 
where 


(3.7) 


(Recall  X(p,q)  *  U)^(p,q)  for  ±q  >  0.)  To  discover  the  behavior  of 
v^(t,x,y)  for  (x,y)  S  R2(d)  and  t  -*■  00  it  will  be  convenient  to  write 
<J>_(x,y,p,q)  in  a  way  that  puts  in  evidence  its  behavior  for  y  -*•  ±°°.  To 
this  end  recall  that  by  [12,  (1 .53)-(l . 60) ] 

(3.8)  <Mx,y,p,q)  =  (2tt) “ 1  c(»)  (2q)  1/2  eip  X  ^(y.  |p|»X) 

for  (p,q)  =  X+(p,X)  €  C+  (and  hence  X  =>  X(p,q)  ■  c2(“)(|p|2  +  q2)). 
Moreover,  by  [12,  (4.5),  (4.6)  and  (8.41)]  one  can  write 


where 


(3.10)  T+(y,p,X)  -  T+(y,X)  4>„(y,U,X)  exp  {iy  q^(y,X)}  -*-T+(y,X)  ,  y  -*■  -®. 

Similarly,  by  [12,  (4.1)  and  (4. 6)-(4.11) ]  one  can  write 


(3.11) 


,  00  ]  1/2 

♦+(y»U.*)  -  fcffi^x)  lI+^y*M,X)  exp  t-iyq+(y,X)} 
+  R+(y»P»X)  exp  (iyq+(U,X)}] 


where 

i+(y»u,X)  -  4>2(y,y,X)  exp  Uyq+(u,X)}  -*■  l 
(3.12)'  •  y  -*■  +». 

R+(y»y»X)  -  R+(U,X)  <My.U,X)  exp  {-iy  q+(y,X)}->-  R+(y,X) 


Combining  (3.6),  (3.8)  and  (3.11)  gives 
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VjU.x.y) 


c(°°)  p(°°)  1/2  ( 2tt) ~ 3/2 


exp  {i(x  •  p+yq  -  tw+(p,q)  }  I+(y,  |p  j ,  X)  £.(p,q)dpdq 


(3.13) 


+  c(°°)  P(°°)  1/2  (2ir)  3/2  j  exp  (i(x  *  p- yq  -  tu+(p,q)  }  R+(y,  |p  | ,  X)  £(p,q)dpdq . 

*  r 


It  is  natural  to  expect  that  in  R+(d)  the  partial  wave  v^(t,x,y)  will 
propagate  as  t  +  ®  into  regions  where  y  is  large  and  hence  I+(y,  |p|,X) 
and  R+(y, |p|,X)  are  near  their  limiting  values.  Thus  the  representation 

(3.13)  suggests  the  conjecture  that 

(3.14)  v^t,*)  ~  +v^(t,«)  inL2(R3(d)),  t 

where  v£  and  v^  are  defined  by 


(3.15)  v£(t,x,y)  =  c(°°) p(°°)  :^2  (2tt)  3^2  j  exp  {i(x*p+y  q  -  tu>+(p,q)  K(p,q)dpdq 


v£(t,x,y)  -  c(°°)p(®)1/2  ( 2tt) “ 3/2  j  exp  (i(x*p-yq  -  tu>+(p,q)}  x 


(3.16) 


x  VlPl>A)  Jl(P»<l)dPdcl 

-  c(»)p(»)  1/2  ( 27T) ~ 3/2  |  exp  {i(x*p+y  q  -  tu)+(p,q)>  x 

-C+ 

x  R+(|p| ,X)  £(p,-q)dpdq 


where  -C+  -  { (p,q)  :  (p,-q)  €  C+}  -  {(p,q)  :  q  <  -a| p| } .  Note  that  v° 
and  v^  are  waves  in  a  homogeneous  medium  with  density  p  (<*>)  and  sound 
speed  c(«°)  .  More  precisely, 
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(3.17) 


Vp(t, *) 


v*(t,-) 


exp  (-i  t  c(“)  Ai/2)  hc 


exp  (-i  t  c(“)  k1,z)  hi 


where  A0  is  the  self adjoint  realization  in  L2(R3)  of 

-A  ■  -(32/3xJ  +  32/3x2  +  32/3y2)  and  h^  and  h^  are  the  functions  in 

L2(R3)  whose  Fourier  transforms  are 


(3.18) 


h^p.q)  *  c(»)  p(“)  1/2  f.(p,q)  -  c(°°)  p(°°)  U2  X^P.q)  h_(p,q), 


1/2 


h^(p,q)  -  c(«)  p(“) 1/2  R+(|pj,X)  2(p,-q) 


c(»)  p(») 1/2  R+( jpj.X)  (1  -  X+(P»q))  h_(p,-q) 


Both  functions  are  in  L2(R3)  because  h_  £  L2(R3)  and  lR+(|p,A)|  <_p1/,2(«>)  by 

A  t 

the  conservation  law  [12,  (1.46)].  Moreover,  supp  h^  C  -C+  and  hence 
the  theory  of  asymptotic  wave  functions  for  d'Alembert's  equation 
[9,  Ch.  2]  implies  that  v^(t,*)  -  0  in  L2(R^(d))  when  t  -*•  00 .  Combining 
this  with  (3.14)  gives 


(3.19)  v^(t,*)  ~  v£(t,*)  in  Lz(R|(d)),  t  ■+  °°. 

Now  consider  the  behavior  of  v^(t,x,y)  for  (x,y)  e  R3(d),  t 
Combining  (3.6),  (3.9)  and  (3.10)  gives 


VjU,x,y)  -  c(«)  p(“)1/2  (2u)-3/2  j  exp  {i(x*p+  yq_  -  tu>+) }  * 


(3.20) 


*  T+(y,  |p  ,A)  £.(p,q)dpdq 


where  q_  -  q_(|p|,A),  w+  -  u)+(p,q)  and  X  -  A(p,q)  -  u>+(p,q).  This 
representation  suggests  that 
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(3.21) 


v5(t,*)  ~  v»(t,*)  in  L2(R3(d)),  t  -*■  “ 


where 


(3.22) 


v|(t,x,y)  -  c(»)  p(°°)1/2  ( 2tt) “ 3/2  j  exp  (i(x  •  p  +  yq_-  tco+)}  * 


r+(  i  P  |  »*)  *-(p,q)dpdq 


Now  the  mapping  (p,q)  ■+•  (p.q1)  =*  X'(p,q)  ■  (p>q_( |p| >w£(p,q)))  with 
domain  C+  has  range  X' (C+)  =  =  R+(0),  Jacobian  3(p,q)/3(p,q ') 

m  c2(-°°)q,/c2(00)q  and  satisfies  w+(p>q)  =  a)_(p,q').  Thus  (3.22)  implies 


the  representation 


Vo(t,x,y)  -  c(°°)  p(°°)l/2  (2ir) - 3/2  f  exp(i(x*  p  +  yq'  -tu>  )}  * 

*  Jr* 


(3.23) 


*  T+(|pj  ,u>2(p,q'))  &(p,q) (c2 (-°°)q  '/c2 (°°)q)dpdq  ' 


where  q  ■  q(|p|,q')  *  /a2 (  p  2  +  q'2)  +  q'*  .  Note  that 


(3.24) 


v|(t,*)  -  exp  (-1 1  c(-®°)  Aj/2)  h£ 


where  h|  e  L2(R3)  has  Fourier  transform 


(3.25)  h£(p,q')  -  c(«)  p(»)  1/2  T+(  |p|  ,U2(p,q'))  Jl(p,q)  (c2(-«)q  7c2  («)q)  . 

Since  supp  h£  C  R^  the  results  of  [9,  Ch.  2]  imply  that  v|(t,*)  -  0  in 
L2(R^(d))  when  t  ■+•  <*>.  Combining  this  with  (3.21)  gives 


(3.26) 


v9(t,*)  -  0  in  L2(R3(d)),  t  -*■  «. 


Analogous  conjectures  concerning  vm(t,*)  and  vn(t,*)  will  now  be 
formulated.  Only  the  main  steps  of  the  calculations  will  be  given  since 
the  method  is  the  same  as  for  v«(t,»). 
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Behavior  of  v  .  v  has  Che  representation 
_ m  m 

(3.27)  vm(t,x,y)  =  4>_(x,y,p,q)  exp  (-1 1  U)+(p,q))  m(p,q)dpdq, 

Co 

by  (3.5),  where 

(3.28)  <j>_(x,y,p,q)  -  (2ir)-lc(“)  (2q)  eip‘x  i|»0(y,  |p|  ,X) 

for  (p,q)  -  X0(p,X)  S  C0  (and  hence  X  *  X(p,q)  ■  u>£(p,q)).  Moreover 
(see  [12,  (4.18)-(4.24)J), 


q/oo1)  1 1/2 

(3.29)  ^(y.U.X)  -  T0(y,u,X)  exp  (yqUy.X)) 


where 


(3.30)  T0(y,y,X)  -  T0(y,X)  <j>3(y,y,X)  exp  (-yq^(y.X))  -*•  T0(y,X),  y  -*■  -» 


and 


(3.31) 


where 


4»0(y,y,X)  -  '  Uo<y»M,A)  exp  {-lyq+(y,x>} 

+  R0(y,y,X)  exp  {iyq+(y,X)>] 


(3.32)  \ 


I„  (y»y*X)  -  4>2(y,y»X)  exp  {iy  q+(y,X)}  -*■  1 


R0 (y.y.X)  -  R„(y,X)  <()j(y,y,X)  exp  {-iy  q+(y,X)}-^  R0(p,X) 
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Combining  (3.27),  (3.28)  and  (3.29)  gives 


vm(t,x,y)  =  c(°°)  pO*’)1'*2  (2tt)  3/2 


exp  (i(x  •  p  -  tw+) }  x 


(3.33) 


xT0(y,|p|,X)  exp  (y  q^)  m(p,q)dpdq 


where  u)+  -  oo+(p,q),  X  *  ui£(p,q)  and  q^  *  q_^ ( | p | » A)  .  Since 


(3.34) 


T  (y,  1  p | , X)  exp  (y  q^C  | p [ , X) )  -*-0,  y 


equation  (3.33)  suggests  that 


(3.35) 


vm(t,*)  -  0  in  L2(R3(d)),  t  -*■  ®. 


Similarly,  combining  (3.27),  (3.28)  and  (3.31)  gives 


vm(t,x,y) 


c(°°)  p(°°)1/2  (2tt) 


»/2  (">^"3/2 


exp  {i(x*p+yq-  tw+)}  x 


x  Io  | P { »^)  m(p,q)dpdq 


(3.36) 


+  c(°°)  p(“) 3/2  (27 r)  3/2  |  exp  {i(x  •  p-yq  -  tu)+)} 


x  R0(y,|p  ,X)  m(p,q)dpdq 


which  suggests  that 


(3.37)  vm(t,«)  ~  v«(t,*)  +  v3(t,*)  in  L2(R].(d)),  t  -  ® 


where  v®  and  v3  are  defined  by 
mm 


(3.38) 


v£(t,*)  -  exp  (-itc(»)  A3/2)  hffl 


v3( t,*)  -  exp  (-itc(»)  Ai/2)  h3 


> 
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and  h  and  hl  are  the  functions  in  L2(R3)  whose  Fourier  transforms  are 
mm 

hm(p.q)  *  c(«)  p(®)l/2  m(p,q)  -  c(»)  p(«>)  1/2  Xo (P»q)  h_(p,q), 

(3.39)  ■ 

h^(p,q)  -  c(»)  p(°°)1/2  R0  ( |  p |  ,X)  m(p,-q) 

-  c (°°)  p(°°)  1/2  R0  (  |  p|  ,X)  (1  -  Xo  (p»q))  h_(p,-q)  . 

Note  that  supp  h1  C  R3  and  hence  v1(t,*)  ~  0  in  L,(R?(d))  when  t 
*  m  —  m  *• 

Combining  this  with  (3.37)  gives 

(3.40)  vm(t,‘)  -  v£(t,  •)  in  L2(R^(d)) ,  t  -*• 

Behavior  of  v  .  v  has  the  representation 
_ n  n 

(3.41)  v  (t,x,y)  *  <{>_(x,y,p,q)  exp  (-itoj  (p,q))  n(p,q)dpdq, 

;c_ 

by  (3.5),  where 

(3.42)  <J>_(x,y,p,q)  -  (2TT)_1c(-«)(2|q|)1/2  eip‘X  K«.(y,|p|,X) 

for  (p,q)  ■  X_(p,X)  €  C_  (and  hence  X  ■  X(p,q)  *  io2(p,q).  Moreover 
(see  [12,  (4.5)  and  (4 .12)-(4 .17) ]) , 

(3.43)  ^_(y,U»X)  -  T_(y,p,X)  exp  Uyq+(y,X)} 

where 

(3.44)  T_(y,ja,X)  ■  T_(y,X)  <J>,  (y,U,X)  exp  {-iy  q+(y ,X) }  -*  T_(u,X) ,  y  -*•  « 
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and 


(3.45) 


^_(y,U,X)  *  [^^_]  [I_(y,y,X)  exp  {iy  q_(y,X)} 

+  R_(y,u,X)  exp  (-ly  q_(y,X)}] 


where 


i_(y.u,X)  -  <J>3(y,y,X)  exp  {-iyq_(y,X)>  -»•  l 


(3.46)  \ 


R_(y.V»^)  -  R_(V,X)  <J>4(y,y,X)  exp  {iy  q_(y,X)}->-R_(y,X) 


Combining  (3.41),  (3.42)  and  (3.45)  gives,  after  simplification  using 
q_(|p|»(*JI(p»q))  -  (q2)12-  -q  for  (p,q)  €  c_. 


vn(t,x,y) 


(3.47) 


c(-®)  p(-“)  l^2  (21T)  3^2  exp  {i(x  •  p+yq  -  tin  ) }  * 

*  C 


*  I_(y»|P  »*)  n(p,q)dpdq 


+  c (~09)  p(-°°) 1,1  ( 2ir) ~ 3/2  exp{i(x*  p-yq  -  tw  )} 

JC_ 

x  R_(y* | P | »X)  n(p,q)dpdq . 


This  suggests  the  asymptotic  behavior 

(3.48)  vQ(t,*)  ~  vq(c»  *)  +  v^(t,  *)  in  L2(R3(d)),  t  ->•  ®, 


where 


(3.49) 


v®(t,*)  -  exp  (-itc(-®)  Aj/2)  hQ 
v^(t,*)  -  exp  (-i  t  c(-»)  Aj^2)  h^ 
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and  h  and  h1  are  the  functions  whose  Fourier  transforms  are 
n  n 


(3.50) 


hn(p,q)  =  c(-“)  p(-°°)  1/2  n(p,q)  =  c(-“)  p(-°°)  1/2  X_ (P> 9)  h_(p,q). 


|  h^(p,q)  *  c(-°°)  p(-°°)  1/2  R_ ( } p |  ,X)  n(p,-q) 

=  c(-“)  p(-oo)  !/2  R_(|p|,X)  (1  -  x_(p»9))  h_(p,-q)  . 

Note  that  supp  h^  C  and  hence  v^(t,*)  ~  0  in  L2(R3(d))  when  t  ■+•  00 . 
Combining  this  with  (3.48)  gives 

(3.51)  vQ(t,  •)  ~  v°(t,  •)  in  L2(R^(d))  ,  t  - 

Finally,  combining  (3.41),  (3.42)  and  (3.43)  gives 


vn(t,x,y) 


=  c(-°°)  p(-“)  1/2  (2tt)_ 3/2 


exp  {i(x  •  p  -  y  q+ -  tai_)  }  x 


(3.52) 


x  T_Cy » 1 P | , X)  n(p,q)dpdq 

where  ■  oa_(p,q),  X  *  u)2(p,q)  and  q_j_  ■  q+(|pj»X).  This  suggests  that 

(3.53)  vn(t,  *)  ~  v2(t,*)  inL2(R|(d)),  t  * 


where 


v2(t,x,y) 

(3.54) 


c(-°°)  p(-°°) 1,2  (  27T)  —  3/2 


exp  (i(x  •  p  -  y  q+  -  tu>_ 


)} 


x  T_( |p| ,X)  n(p,q)dpdq. 


Now  the  mapping  (p,q)  (p,q')  -  X"(p,q)  -  (p,-q+(p,X(p,q)))  maps  C_ 

onto  X"(C_)  ■  -C+,  has  Jacobian  3(p,q)/3(p,q ')  ■  c2  (°°)q '/c2 (-°°)q  and 
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satisfies  w_(p,q)  *u>+(p,q').  Thus 


(3.55) 


v*(t,‘)  *  exp  (-i  t  c(°°)  Aj/2)  h2 


where  h2  has  the  Fourier  transform 
n 


h*(p,q')  *  c(-“)  P(-°°)  1/2  T_(|p[,w£(p,q'))  n(p,q(p,q')) 


(3.56) 


x  (c2  (°°)q  '/ c2(-°°)q) 


Moreover,  supp  h*  C  R3  and  hence  v^(t,*)  ~  0  in  L2(R^(d)),  t 
Combining  this  with  (3.53)  gives 


(3.57) 


v  (t,*)  ~  0  in  L2(R^(d))t  t  -►  “. 


The  asymptotic  behavior  of  v^(t,*)  for  t  -*•  00  may  be  obtained  from 
the  three  cases  analyzed  above  by  superposition,  equation  (3.4).  Thus 
equations  (3.19),  (3.26),  (3.35),  (3.40),  (3.51)  and  (3.57)  imply 


vj[(t,‘)  +  v“(t,')  in  L2(R|(d)) 


(3.58) 


vf(t,*) 


vj(t,-) 


in  L2(Rf(d)) 


t  0°. 


On  combining  this  with  the  definitions  of  v£,  v£  and  v£,  equations  (3.17), 
(3.18),  (3.38),  (3.39),  (3.49)  and  (3.50),  one  is  led  to  formulate 
Theorem  3.1.  For  every  h  6  5C  let  v£(t,»)  be  defined  by 


(3.59)  v°(t,x,y) 


[exp  (-itc(<*>)  Aj/2)  h+(x,y),  (x,y)  €  R^.(d) , 


exp  (-itcH  A,/2)  h“(x,y),  (x,y)  e  Rf  (d) , 


where  h  and  h~  are  the  functions  in  L2(R3)  whose  Fourier  transforms  are 


given  by 


1 
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(3.60) 


h+(p,q) 


c(oo)  p(«)  1/2  h_(p,q)  ,  (p,q)  €  R_J., 


,  (p, q)  €  R3, 


and 


(3.61) 


_ 

h  (p,q) 


0  ,  (p, q)  £  R+, 

c(_o°)  p(-°°)  1/2  h_(p,q)  ,  (p , q)  SR3. 


Then 


(3.62) 


lim  8 vf (t, •)  -  v°  (t,  •)& 

£-»CO 


X 


0. 


Theorem  3.1  implies  corresponding  asymptotic  estimates  for  the 
free  component  uf(t,*)  =  Pfu(t,‘)  =  Re  {vf(t,*)}  of  the  acoustic  potential 
u(t,*).  Indeed,  if  u^(t,*)  is  defined  by 


(3.63)  u°(t,*)  »  Re  (vjUt,*)} 

then  Theorem  3.1  and  the  elementary  inequality  |Re  z|  <  |zj  imply 
Corollary  3.2.  For  all  h  e  X  one  has 


(3.64) 


lim  i  u, (t,  •) 


(t,  =  0. 


If  the  initial  state  h  has  derivatives  in  Jf  then  u^(t,*)  and 
u°(t,*)  have  the  same  derivatives  in  X  and  (3.64)  can  be  strengthened  to 
include  these  derivatives.  In  particular t  one  has 

Corollary  3.3.  For  all  h  €  L*(R3)  =■  D(A1/2)  one  has 


(3.65) 


lim  ||D.  u. ( t ,  • )  -  D.  u^t,*)!^  =  0,  j  =  0,1, 2, 3, 
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where  D0  *  3/3t,  Dj  “  3/3xj,  D2  *  3/3x2  and  D3  =  3/3y.  Equation  (3.65) 
is  equivalent  to  convergence  in  energy: 

(3.66)  lim  E(u^  -  u£,R3,t)  =  0. 

t-K» 

Corollary  3.3  can  be  proved  by  applying  the  method  of  this 
section  to  the  derivatives  v^(t,*)  (j  ■  0,1, 2, 3)  which  are  given  by 

integrals  of  the  same  form  as  (3.2).  Detailed  proofs  for  the  case  of  the 
Pekeris  profile  were  given  in  [10]. 

Proof  of  Theorem  3.1.  The  remainder  of  this  section  is  devoted 
to  the  proof  of  Theorem  3.1.  The  decomposition  (3.3)  is  used  for  the 
proof.  Moreover,  for  brevity,  only  the  asymptotic  equality  (3.53)  for 
vn(t,*)  is  proved.  The  remaining  five  cases,  namely  (3.14),  (3.21), 
(3.35),  (3.40)  and  (3.48)  can  be  proved  by  the  method  used  for  (3.53). 

As  a  first  step,  (3.53)  will  be  proved  for  the  special  case  of  n(p,q) 

€  C0(CJ,  the  set  of  continuous  function  with  compact  supports  in  the 
open  cone  C_.  The  general  case  will  then  be  proved  by  using  the  fact 
that  C0(C_)  is  dense  in  L2(C_). 

For  functions  n(p,q)  €  C9(C_)  the  integrals  defining  v^  and  v* 
converge  point-wise,  as  well  as  in  JC,  and  one  can  write 


where 
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(3.68)  w(y,p,t) 


■o 

exp  (-i(yq+  +  tu>_)}  [T_(y,  |  p  { -  T_(  |p|,caf)  )n(p,q)dq 

— 00 


Parseval's  formula  in  L2(R2),  applied  to  (3.67),  gives 


(3.69) 


vn(t,x,y)  -  v2(t,x,y) j2dx 


c2(_<»)  p(-oe)(27r)~l  (  |w(y,p,t)  |2  dp. 

jR2 


On  integrating  this  over  y  >  d  one  finds 

(3.70) 

The  last  relation  implies  that  to  prove  (3.53)  it  is  sufficient 
to  prove  that 

(3.71)  w(-,t)  -  0  in  L2(R|(d)),  t  -  ». 

To  this  end  it  will  be  convenient  to  change  the  variable  of  integration 
in  (3.68)  from  q  to  w  ■  w_(p,q)  ■  c(-°°)  /jp  | 2  +  q2 .  Solving  this 
equation  for  q  <  0  gives  q  *»  -(u)2  c  2(-°°)  -  |p|2)^2  *  ~q_(|p|,u>2)  with 
oj  >  c(-°°)|pj.  Hence  (3.68)  can  be  written 


(3.72) 


where 


w(y,p,t) 


IpI 


exp  (-itu>)  W(y,p,u>)du> 


W(y,p,u)) 

(3.73) 


-2 

C 


(-00)  exp  {-iy  q^.(  |p  |  ,<*>2)  }[T_(y,  ]p  |  .uj*)  -  T_(  jp  |  ,uiJ)  ]n(p,-q_(  |p|  ,u)2))w 


q_(  jpj  ,cj2) 
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> 


I 


I 


1 


♦ 


» 


» 


► 


The  assumption  that  n  e  C0(C_),  together  with  (3.44),  implies  that 
W  £  C0(R  x  F)  where  T  »  {(p,w)  :  <u  >  c(-°°)|pj}.  Moreover,  by  a  standard 
partition  of  unity  argument,  one  may  assume  without  loss  of  generality 
that 

(3.74)  supp  W(y,*)  C  {(p,w)  :  |p|  <  p0  and  0  <  w0  <  u>  <  Wj } 
for  all  y  €  R  where  ojo  >  c(-«)pa.  This  in  turn  implies  that 

(3.75)  w(y,p,t)  -  I  1  exp  (-itiu)  W(y,p,u)daj 

H 

and 

(3.76)  supp  w(y, • , t)  C  B(p0)  -  (p  :  |p|  <  pQ} 
for  all  y  €  R  and  t  €  R.  Thus 


(3.77) 


lv(-,t)|f  /R3/d-.x  -  I  f  |w(y,p,t)  ( 2  dpdy 

*d  JB(p0) 

-  [  [  |w(y,p,t) |2  dpdy 

Jd  JB(p0) 


+  J  [  |w(y,p,t)|2  dpdy 

Jy0  yB(p») 


for  any  yfl  >  d.  The  proof  of  (3.71)  will  be  derived  from  (3.77)  and  the 
following  two  lemmas. 

Lemma  3.4.  Let  n€  C0(C_)  and  assume  that  (3.74)  holds.  Then 
for  each  d  €  R,  y0  >  d  and  p0  >  0  one  has 
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(3.78)  lim  w(y,p,t)  -  0, 

t-+co 

uniformly  for  all  (y,p)  S  [d , y 0 ]  x  B(p0)  . 

Lemma  3.5.  Under  the  hypotheses  of  Lemma  3.4  there  is  a  constant 
C  *  C(n)  such  that 

(3.79)  |w(y,p,t)|  <  C  y 1_0t 

for  all  y  >  0,  p  £  B(p0)  and  t  e  R,  where  a  >  3/2  is  the  constant  of 
condition  (1.6). 

Proof  of  Lemma  3.4.  The  proof  is  based  on  a  well-known  proof  of 
the  Riemann-Lebesgue  lemma.  Note  that  by  (3.75)  one  has 


w(y,p,t) 


(3.80) 


■a)1-(7r/ 1) 

exp  (-iw  t)  W(y,p,aH-(Tr/t))du 

^u0-(ir/t) 


exp  (-i  t  w)  [W(y,p,ci))  -  W(y,p,w+(iT/t))  ]dco 


1  f^0 

exp  (-itaj)  W(y,p,oH-(ir/ t))do> 

^u)0~  (it/  t) 


+  4  j  exp  (-i  t  m)  W(y,p,w+(iT/t))dw. 

The  limit  relation  (3.78)  is  obvious  from  (3.80)  and  the  continuity  of 
W.  The  uniformity  of  the  limit  follows  from  (3.80)  and  the  uniform 
continuity  of  W  on  compact  subsets  of  R  *  T. 

Proof  of  Lemma  3.5.  Note  that  by  (3.72),  (3.73)  and  (3.44)  one 


has  the  estimate 


I 


9 


9 
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(3.81) 


|w(y,p,t) 


(y>  |p|.<*>2)  exp  {-iy  q+(|p|,w2)}-  l|  M(p,u))dw 


for  all  y,  p  and  t  where 


(3.82) 


M(p,ai)  ■  c“2(-°°)  |T_(|pj,U)2)  | 


n(p,-q_(  |p  |,0)2))cj 


q_(  IpI.w2) 


is  continuous  for  |p|  <  p0,  w0  <  to  <  toj.  It  follows  that  for  any  y  €  R, 

P  e  B(p0)  and  t  e  R 

(3.83)  |w(y,p,t)  |  <  M„  sup  |$i  (y,  |p|  ,to2)  exp  {-iy  q+(  |p  |  ,to2)  }  -  l| 

where  M0  ■  M„(n)  *  sup  M(p,to)  and  the  suprema  are  taken  over  all  |p(  <  p0 
a«.d  u)0  <  to  <  Uj.  The  proof  of  (3.79)  will  be  based  on  (3.83),  conditions 
(1.3),  (1.4),  (1.6)  on  p(y)  and  c(y)  and  the  proof  in  [12]  of  Theorem 
2.1;  see  [12,  p.  27ff].  Note  that  because  of  the  continuity  of  w  it  will 
suffice  to  prove  (3.79)  for  all  y  >  yt  where  -  yj(n)  is  a  positive 
constant. 

The  solution  ^(y.y.A)  with  X  >  c2(-°°)y2  >  c2(°°)y2  satisfies 
[12,  p.  27ff ] 

(3.84)  4>i(y,M,X)  -  exp  (iyq+(y,A)}  (nx  +  n2) 

where  n  -  (n1»n2)  is  characterized  on  y  >  yx  as  the  unique  solution  of 
the  integral  equation  [12,  (2.64)]  which  can  be  written 


(3.85) 


n  -  n®  +  K(y,X)n,  n®  -  (i,o). 


The  kernel  K(y,y',y,A)  -  (K^ (y,y* ,y,X))  is  defined  by 


» 
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(3.86) 


and 


KjjCy.y'.y.*) 


0  ,  yx<y'  <  y, 


^^(y'.PfX),  y'  >  y. 


(3.87)  K  (y.y'.p.A)  = 


o  »  yx  <  y'  <  y. 

|-exp  {-2i (y-y ')<!_!_ (u,X)  E^  (y' ,y,X) ,  y'  >  y. 


where  [12,  p.  27] 


(3.88) 


E(y,y,A)  -  B-1(y,A)  N(y,y, A)  B(y,A). 


From  these  relations  one  has  [12,  (2.64)] 


<J>i(y.y,A)  exp  (-iyq+(M,X)}  -  1  ■  Hj  +  U2  -  1 


(3.89) 


#<»  j-OO 

-]  Eij(y,,y,X)  n j (y ' ) dy ’  -  exp  {-2i(y-y’)q+(y,A)}  E^  (y ’ ,y,A)n.j  (y ’)dy ’ 


and  hence 


2  *00 

(3.90)  |<j>1(.v,y,X)  exp  {-iy  q^y.A)}- l|  <  f  j  |E.k(y’ ,y,X)  |  |nk(y')  |dy  ’ . 

j,k-l  }y 

Using  [12,  (2.66)]  and  the  continuity  of  B(y,X)  on  X  >  c2(-»)y2  it  can 
be  shown  that 

(3.91)  DK(y,A)l  <  1/2  for  0  <  y  <  p0 ,  u>2  <  X  <  uj2 
provided  yt  -  yt(n)  is  large  enough.  Thus 
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(3.92) 


!ns  < 


-  1  -  IKI  1  -  IKI  - 


<  2 


for  0  <  y  <  p0,  uij  <  X  <  wf .  Combining  this  with  (3.90)  and  [12,  (2.66)] 
gives 


(3.93) 


((My.U.A)  exp  {-i y q, (p,X)}  -  l|  <2  \  |Ek(y',y,A)j  dy' 

<  C,  j  | p(y * )  -  p(“) |  dy' 

sy 

+  c2  J  |c(y ')  -  c(°»)  I  dy' 

*  XT 


for  the  same  values  of  y  and  A  where  Cx  and  C2  depends  only  on  p0,  w0 
and  ial.  It  follows  from  (3.93)  and  (1.6)  that 

(3.94)  sup  |<j>j  (y»  |p|  ,w2)  exp  {-iy  q+(  |p|  ,w2)>  -  l|  <  C3  y1_Ct 

for  all  y  >  yx(n)  where  C3  depends  only  on  p„,  u>„  and  (i.e.,  n) . 
Combining  (3.94)  with  (3.81)  gives  (3.75). 

Proof  of  Theorem  3.1  (completed) .  Lemma  3.5  implies  that 

(3.95)  f°  [  |w(y,p,t)  j 2  dpdy  <  (irp§  Cz/2a  -  3)  y\  2a 

y0  JB(Po) 

where  3  -  2a  <  0.  Thus  given  any  e  >  0  there  is  a  y0  m  y0(e,n)  such  that 


(3.96) 


B(P0) 


|w(y»p»t)  | 2 


dpdy  <  e 


for  all  t  €  R.  Equation  (3.77),  together  with  Lemma  3.4  and  the 
estimate  (3.96),  imply  that 
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(3.97) 


lim  sup  II  w(* ,  t)  II  2  <  e. 


L2(R*(d))  - 


Since  e  >  0  is  arbitrary  this  implies  (3.71),  as  required. 


The  arguments  given  above,  applied  to  v^,  vm  and  vn,  show  that 
the  conclusions  of  Theorem  3.1  hold  for  all  h  such  that 
h_  €  C0(C+  u  C0  u  C_) .  Moreover,  this  set  is  dense  in  L2(R3)  and  hence 
C0(C+  UC0U  C _)  is  dense  in  Xf  =>  Pf  X  by  [12,  Cor.  8.12].  These 
facts  can  be  used  to  extend  (3.62)  to  all  he  X  because  the  mappings 
U(t)  :  Xf  f  L2(R3)  and  U„(t)  :  3Cf  -*■  L2(R3)  defined  by  U(t)hf 
*  exp  (-It  A1^2)  h^  and  U0(t)h^  =  v^(t,*)  are  uniformly  bounded  for  all 
t  €  R  (see  [11,  p.  32]).  The  density  argument  needed  to  extend  (3.62) 
to  all  h  €  JC  has  been  given  in  many  places;  see,  for  example,  [9,  Ch.  2] 
or  [10,  p.  260]. 
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§4.  Transient  Guided  Waves. 

The  asymptotic  behavior  for  t  -*■  00  of  the  guided  component 
Ug(t,*)  ■  Pg  u(t, *)  is  derived  in  this  section.  u^(t,*)  is  a  sum  in 
of  mutually  orthogonal  partial  waves  u^Ct,*)  ■  u(t,*)  *  Re  (v^(t,*)}, 

1  <  k  <  N0.  The  starting  point  for  the  analysis  is  the  integral 
representation 


(4.1) 


vk(t,x,y)  -  |  ^k(x,y,p)  exp  {-it wk(|p|)  hk(p)dp 


where 


(4.2) 


hk(p)  ■  |  3  \(x»y»P)  h(x,y)  c-2 (y)  p-1(y)dxdy 


and  the  integrals  converge  in  30  and  Lz(^k),  respectively.  The  integral 
in  (4.1)  can  be  written 


(*•3)  vk(t,x,y)  »  ~  |  exp  {i(x  •  p -  t  oik(  jp | ) ) }  </>k(y,p)  hk(p)dp. 


This  is  an  oscillatory  integral  that  can  be  estimated  by  the  method  of 
stationary  phase  when  6  C*(flk) .  To  apply  the  method  define 


r  »  /t2  +  |  x  I 


(4.4) 


t  -  r  $o.  *i  ■  r  x2  -  r  52 


C  -  (e0,Cx,Ca)  €  S*  C  R3 
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where  S2  denotes  the  unit  sphere  in  R3.  Then  (4.3)  takes  the  form 


(4.5) 


vfc(t,x,y) 


exp  {ir  0k(p,C)}  gk(p»y)  dp 


where 


(4.6) 


6k(p.S)  -  ZiPi  +  C2Pa  -  S0  uk(lpl) 
gk(p»y)  *  (2ir)-1  i|>k(y,p)  h^p) 


Estimates  for  large  r  of  the  integral  in  (4.5)  are  needed  that  are 
uniform  for  (5,y)  in  compact  subsets  of  S2  *  R.  Such  estimates  are 
provided  by  a  version  of  the  method  of  stationary  phase  due  to 
M.  Matsumura  [5].  A  form  of  Matsumura's  results  applicable  to  (4.5)  was 
presented  in  [11,  Appendix],  This  result  is  applied  below  to  estimating 

V 

The  phase  function  9k  has  a  point  of  stationary  phase  if  and  only 
if 

(4.7)  Uk(|p|)p/|p|  -  x/t 


where 


(4.8)  Uk(!p|)  -  u)^(|p|) 

is  the  group  speed  associated  with  the  dispersion  relation  u  -  cuk(|p|). 
It  will  be  assumed,  for  brevity,  that  Uk(p)  is  a  monotone  decreasing 
function  that  maps  0k  onto  (cm,c(so)).  In  this  case  (4.7)  has  a  unique 
solution  if  |x|/t  lies  in  the  range  of  Uk(|p|);  that  is 
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(4.9) 


cm  <  |x|/t  <  c(°°). 


and  no  solution  otherwise.  The  solution  is  given  by 


(4.10) 


P  ■  Qv( lxl/t)  x/ ix 


where  is  the  inverse  function  to  U^.  By  calculating  the  Hessian  0£ 


one  can  show  that 


(4.11) 


|det  0«p,C)|  *  t2  U.  (jp|)  |  U*  (  t  p  | )  |/|p| 


and  sgn  9k(p,5)  ”  0.  In  particular,  each  point  of  stationary  phase  is 
non-degenerate  and  makes  a  contribution 


(4.12)  v.  (t,x,y,p)  - 


(p|J/2exp  {i(|xj  |p|-to»k(|p|))}  i^k(y,|p|)hk(p) 

(|p|)|u'(|P|)|}1/2 


to  the  integral  in  (4.3),  where  p  is  given  by  (4.10).  For  |x|/t  outside 
the  interval  (4.9)  there  is  no  point  of  stationary  phase.  Thus  the 
stationary  phase  approximation  to  vk(t,x,y)  is  given  by 


(4.13) 


vk(t,x,y)  •  x(lx|/t)  v"(t,x,y,Qk( jx|/t)  x/|x|) 


where  x  is  the  characteristic  function  of  the  interval  (cm,c(°°))  and  one 


Theorem  4.1.  For  all  h  e  3f  such  that  h^  6  C(J(fik)  there  exists  a 
constant  C  *  Ck(h)  such  that 


(4.14) 


|vk(t,x,y)  -  vk(t,x,y)|  <  C/tJ 


for  all  t  >  0,  x  £  R2  -  (0)  and  y£  R. 


Theorem  4.1  can  be  proved  by  application  of  Theorems  A.l  and  A. 2 


of  [11].  The  proof  that  ^k(y,p)  has  the  required  p-derivatives  is 
lengthy  but  straightforward  and  will  not  be  given  here.  If  h^  is  not  a 
smooth  function  then  Theorems  A.l  and  A. 2  are  not  applicable  and  the 
estimate  (4.14)  may  fail.  However,  the  definitions  (4.12)  and  (4.13)  are 
meaningful  for  all  h  6  X  and  one  has 

Theorem  4.2.  For  all  h  £  J C,  all  t  >  0  and  k  ■  1,  2,  3,*«*  one  has 

(4.15)  v"(t,*)<S3C 

and 


(4.16)  "  lPkh,3C  * 

00 

Moreover,  the  mapping  t  -*•  vk(t,*)  is  continuous  from  R+  to  K  and 

(4.17)  lim  lvk(t,*)  -  v£(t,  •)*£.-  0. 

The  proofs  of  these  properties  are  the  same  as  those  for  the 
Pekeris  profile,  given  in  [10],  and  are  not  reproduced  here.  On 
defining 

(4.18)  u£(t,x,y)  -  Re  (vk(t,x,y)} 


one  also  has 

Corollary  4.3.  For  all  h  €  Jf  and  k-1,  2,  3, • •*, 


(4.19) 


lim  luk(t,»)  -  uk(t,OI3f 
t-*° 


0. 


If  h  €  Lj(R3)  then  uk(t,*)  €  Lj(R3)  and  asymptotic  wave  functions 
for  the  first  derivatives  of  u^  can  be  constructed.  Indeed,  if 
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€  C0(^)  then  the  first  derivatives  of  v^  are  given  by 


Dt  vk(t>x,y) 
(4.20) 

Dj  vk(t,x,y) 
Dy  vk(t,x,y) 


_1_ 

2tt 


2ir 


JL_ 

2tt 


exp  (i(x  •  p-  t  a)k(|p|))}(-iwk(  |p  |))'('k(y,p)hk(p)dp, 


exp  {i(x  •  p  -  t  u^dpl))  (ip..)  ^k(y,p)  hk(p)dp  (j  = 


exp  {i(x  *  p-  t  u^Clpj))}  Dy  ^k(y,p)  hk(p)dp. 


1,2), 


These  integrals  have  the  same  form  as  the  integral  (4.3)  for  vk>  The 
corresponding  asymptotic  wave  functions  are  defined  by 


(4.21) 


vk0(t,x’y’p)  *  (-i  “fcOpI))  vk(t,x,y,p). 


vkj (t,x,y,p)  *  (ipj)  vk(t,x,y,p)  (j  =  1,2), 


vk3(t,x,y,p)  =  Dy  v"(t,x,y,p),  and 

vkj (t»x»y)  *  X( i x|/t)  v” (t,x,y,Qk( |x|/t)  x/|x|) 

for  j  -  0,1, 2, 3.  The  analogue  of  Theorem  4.2  is 

Theorem  4.4.  For  all  h  e  L*(R3),  all  t  >  0  and  k  *»  1,2,3, •••  one 

has 


(4.22) 


vkj(t,-)  6  MR’),  j  -  0,1, 2, 3, 


(4.23)  lv“  <t,-)l*  +  ^  1''kJ<t-->^2(R>,p-d*dy)  '  2 
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and 


(4.24) 


11m  ID 
t-»o 


Vk(t**>  "  Vkj(t*,)BL, 


(R3) 


0,  j  =  0,1,2, 3. 


The  proof  of  Theorem  4.4  is  the  same  as  for  the  special  case  of 
the  Pekeris  profile  which  was  treated  in  detail  in  [10]. 

The  preceding  discussion  was  restricted  to  the  special  case 
where  U^(y)  is  monotonic.  If  U^U)  has  a  finite  number  of  maxima  and 
minima  there  are  a  corresponding  number  of  points  of  stationary  phase 
and  the  form  of  the  asymptotic  wave  function  is  more  complicated  but 
still  tractable.  In  the  case  of  the  Pekeris  profile,  treated  in  [10], 
there  are  two  points  of  stationary  phase.  Cases  that  lead  to 
infinitely  many  stationary  points  have  not  yet  been  encountered.  They 
would  require  additional  analysis. 


§5.  Asymptotic  Distributions  of  Energy  for  Large  Times. 

The  total  energy  of  the  acoustic  field  u(t,x,y),  given  by  (2.14), 
is  constant  for  t  >  T.  The  same  is  true  of  the  partial  waves  u^,  ug  and 
u^,  k  »  1,  2,  3,***.  Moreover,  it  was  shown  in  [12]  that  {P^,Pj ,P2, • ••} 
is  a  complete  family  of  orthogonal  projections  in  X  that  reduces  A.  It 
follows  that 


No-1 

(5.1)  IA1/2  hl£  -  1A1/2  hflj  +  I  |A1/2  ly2. 

k-1 


which  may  be  interpreted  as  an  energy  partition  theorem.  The  partial 
energies 


(5.2) 


E (uf , R3,  t)  =  IA^2  hfl  2, 

E(uk,R3,t)  -  I  A172  t^l2.,  k-1,  2,  3,'“,  NQ-1, 


can  be  calculated  from  the  source  function  f(t,x,y)  and  the  normal  mode 
functions.  The  relationship  between  h  and  f  is  given  by  (2.7),  which 
implies  that 


(5.3)  h_(p,q)  «  i  X-1/2  (p,q)  f_(X1/2  (p,q),p,q) 


where 


(5.4)  f_(w,p,q) 


exp  (iwr)  <P_(x,y,p,q)  f  (t,x,y)  c"2  (y)  p_1  (y)dydxdT 


and  similarly 
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(5.5) 


where 


\(p>  ■  i  \l(|p|)  ^k.(“k(  IpI)»P> 


(5.6)  fk(w,p)  -  ll  s  exp  (iuT)  ^k(x,y,p)  f(T,x, 


y)  c  *(y)  P~ 1 (y)dydxdt . 


It  follows  chat  the  partial  energies  are  given  by 


(5.7)  £(uf,R3.t)  -  (  |f  (X1/2  (p,q),p,q)|2  dpdq 

4  Jb* 


and 


(5.8) 


E(uk,R3,t) 


I  |fk(Vlpl)’P)l2  dp’  k  *  L»  2> 


for  every  t  >  T. 

The  theorems  of  §3  and  §4  make  it  possible  to  calculate  asymptotic 
distributions  of  energy  in  bounded  and  unbounded  subsets  of  R3.  Only  the 
principal  results  are  formulated  here.  The  proofs  are  omitted  since 
they  are  the  same  as  those  for  the  Pekeris  profile  which  were  given  in  [10] . 

The  notation 


(5.9)  Ea0(u,K)  ■  lim  E(u,K,t) 

t-K«> 

will  be  used  whenever  the  limit  exists.  A  first  result  is  the  transiency 
of  all  wavers  with  finite  energy  in  stratified  fluids: 

(5.10)  E*(u,K)  -  0  for  all  compact  sets  K  C  R3. 
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The  Free  Component.  The  results  of  §3  imply  that  in  each  half¬ 
space  R^(d)  (resp.  R*(d))  behaves  like  a  wave  in  a  homogeneous 

medium.  Thus  if  C+  (resp.  C  )  denotes  a  cone  in  R^(d)  (resp.  R^(d)) 
then  Corollary  3.3  and  the  results  of  [9]  imply 


(5.11) 


^(u^.C*)  -  c2(±“)  p(±«)  f  (  |p  1 2  +  q2)  |h_(p,q)  | 2  dpdq 

■  P(±°°)  +  |f_(w±(p,q)  ,p,q)  |2  dpdq. 

J  c~ 


It  follows  that  if 

(5.12)  S  -  {(x,y)  :  dj  <  y  <  d2} 
is  a  slab  then 

(5.13)  E°°  (uf  ,S)  =  0. 

The  Guided  Component.  Consider  the  family  of  cones  defined  by 

(5 .14)  C(e,d)  -  {(x,y)  :  |y  —  d |  <  e|x|} 
where  e  >  0.  Then,  in  contrast  to  (5.11),  one  has 

(5.15)  E°°  (u^.CCe.d))  -  E"(uk,R3)  -  I A 1/2  1^1  J 

for  every  e  >  0  and  k  ■  1,  2,  3,*»*  (see  [10,  Theorem  5.5]).  Finally, 
if  S  is  the  slab  defined  by  (5.12)  then  one  can  show  by  means  of  Theorem 
4.4  that 


J  ^(y.p)  C~2W  P_I(y)dy 


dp. 


(5.16) 
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§6.  Seal-Infinite  and  Finite  Layers. 

The  preceding  analysis  is  extended  In  this  section  to  the  cases 
of  semi-infinite  and  finite  layers  of  stratified  fluid.  The  extensions 
are  based  on  the  normal  mode  expansions  for  these  cases  that  were  derived 
in  [12,  §9].  Only  the  principal  concepts  and  results  are  formulated 
here  since  the  proofs  are  entirely  analogous  to  those  of  the  preceding 
sections. 

Semi-Infinite  Layers.  As  in  [12,  §9]  the  fluid  is  assumed  to 
occupy  the  domain  and  to  satisfy  the  Dirichlet  or  Neumann  boundary 
condition.  Here  the  functions  p(y)  and  c(y)  are  assumed  to  be 
Lebesgue  measurable  and  satisfy 

(6.1)  0  <  pm  <  p(y)  <  PM  <  °°,  0  <  cm  <  c(y)  <  <  00 

and 

(6.2)  |p(y)  ~  P(°°)|  <  C  y^JcCy)  -  c(“)|  <  C  y-® 

for  all  y  >  0  where  Pn,  p^,  p(°°),  cm,  c^,  c(°°),  C  and  a  are  constants  and 

(6.3)  a  >  3/2. 

As  in  [12,  §9]  the  acoustic  propagators  for  p(y)  and  c(y)  corresponding 
to  the  Dirichlet  and  Neumann  conditions  will  be  denoted  by  A°  and  A1, 
respectively.  They  are  selfadjoint  non-negative  linear  operators  in 
*+  -  M*+.c_2(y)  p”1  (y)dxdy)  . 

The  normal  mode  functions  i^(x,y,p,A)  for  A^,  as  defined  in 
[12,  §9],  are  parameterized  by  (p,A)  e  fl  -  ( (p» A)  |  A  >  c2(®)jp|2}. 
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Their  asymptotic  form  for  y  -►  00  is  given  by 


(6.4) 


**< 


x.y.p.q) 


Ip  1  /ai(p*x-yq)  +  Rj  eKp*x+yq) 


2tt 


where  q  -  q(|p|,X)  -  (X  c"2(°°)  -  |p|2)1/5  ,  c(|p|,X)  -  (p(®)/4irq(  |p|  ,X)) 1/2  , 
R^  ■  RJ(|p|,X)  and  (R^(|p|»X)|  -  1.  As  in  the  preceding  sections,  it 
will  be  convenient  to  introduce  new  parameters:  (p,X)  -*■  (p,q) 

■  (p*q(|p|»X))  €  R^  and  normal  mode  functions 

(6.5)  <J>+(x,y,p,q)  -  (2q)1/2c(«)  (|/^(x,y,p,X) 

where 


(6.6)  X  -  X(p,q)  -  c2(")(jp|2  +  q2). 

The  asymptotic  form  of  is 


(6.7)  <(^(x,y,p,q) 


c(°°)p  ^  (°°)  f  i(p.x,qy)  +  Rj  gi(p*x+qy)\ 
(2ir)  Vi  \  J 


y  ^  oo. 


The  second  family  <t»£  defined  by 


(6.8)  <^.(x,y,p,q)  -  4>+(x,y,-p,q) 


is  also  needed.  It  satisfies 


(6.9)  $_(x,y,p,q)  ~ 


c(“)p1/2  («) 
(2tt)  W 


jei(p*x+qy)  +  Rj  gKp'x-qy)^  y 


The  expansion  theorem  of  [12,  §9]  implies  that  the  limits 

-M  f  - 

(6.10)  f|(p,q)  -  L2(R?)-lio  $+(x,y,p,q)  f(x,y)  c“2(y)  p-1(y)dxdy 

1  |x|<M 


tt*oo 


exist.  Moreover,  if 
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(6.16)  v-(t,x,y)  -  <^(x,y,p,q)  exp  (-It  w(p,q))  h^(p,q)dpdq 

1  J  n  3 


where  w(p,q)  *  c(°°)t/|p  | 2  +  q2 .  Moreover,  as  in  §3  one  can  write 


(6.17)  4*^(x,y,p,q) 


c(“)pl/2  (») 


(2ir) 


3/2 


|ei(p-»Ky)  ^(y.p,,)  + 


where 


(6.18) 


lim  1^ (y,p,q)  -  1, 

y^o 


lim  R^(y,p,q)  -  R^(p,X). 

Y** 

Then,  proceeding  as  in  §3,  one  can  prove  the  following  analogue  of 
Theorem  3.1. 

Theorem  6.1.  For  every  h  6  3C+  let  v®(t,’)  be  defined  by 
(6.19)  Vj(t.x.y)  -  exp  (-itc(“)  Aj/2)  h„(x,y),  (x,y)  €  Rj. 

where  h0  €  L2 (R3)  is  the  function  whose  Fourier  transform  is 


(6.20) 


Then 


b0(P»q)  -  1 


c(oo)  p^2  (»)  hJ(p,q),  (p ,q)  €  R|, 
0  ,  (p»q)  €  R^. 


(6.21) 


lim  |vf(t,0  -  v|(t,»)|„  -  0. 


^(y>p,q) 
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Transient  Guided  Waves.  For  both  semi-infinite  and  finite 
layers  the  form  of  the  guided  components  v^(t, •)  is  precisely  the  same 
as  for  the  case  of  an  infinite  layer.  Thus  the  analysis  of  §4  applies 
unchanged  to  these  cases. 
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